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MOLLIER I,S DIAGRAMS FOR COMBUSTION GASES 
IN DATA PROCESSING 

Dr.-Ing. Friedemann Zacharias 


In order to have all the thermal and caloric 
states of combustion gases accessible in a computer, 
closed mathematical approximation equations were 
established for the real factors, the enthalpy and 
the entropy of a real combustion gas. The equations 
approximate the various effects of molecular forces - 
real gas influence and dissociation - at temperatures 
of 200 - 6,000 K, pressures of 0.001 - 1,000 bar, 

and in the range from stoichiometric composition to 
air. Finally, a system of subprograms is listed 
in FORTRAN, by means of which thermodynamic 
calculations can be carried out in the same manner 
as with Mollier I,S diagrams. 


After the publication of the first I,S diagrams for 
gases in 1932 [1], Prof. Pflaum of the Technical University 
in Berlin made the thermodynamic fundamentals of combustion 
engines his objective. In 1960, the second edition of the 
Mollier I,S diagrams for combustion gases [2] was also 
published at his institute. This edition included the most 
modern physical data for the individual gases included in 
the combustion gases. 

With the expanded use of programmed computer systems 
in engineering work, it had become desirable to supplement 
the graphic computing aid of the I,S diagrams for use with 
electronic data processing. 

In the years 1963 to 1965, research work was carried 
out with the objective of representing, by means of equations, 
the thermodynamic functions, which also formed the basis of 
the I,S diagrams. The enthalpy, the entropy, ana the 


1 


real factors of the thermal equation of state for combustion 
gases were developed as closed formulas valid over the entire 
range; 

Temperature: 200 - 6,000 K 

Pressure: 0.001 - 1,000 kp/cm^ (absolute) 

Air ratio: stoichiometric gas - air. 

In a further step, a programming system was established for 
computers to allow them to operate with these formulas in 
the same manner as manual operations with I,S diagrams. 

This decisively simplified, or even made possible for 
the first time, computerized thermodynamic determinations 
on combustion engines, as well as stepwise simulations of 
real processes or optimization computations. 

The Mollier I,S diagrams, though, remain essential 
for a survey of the establishment of thermodynamic computer 
procedures and for quick calculation of individual state 
points and simpler processes. 

The work arranged as follows ; 

1. General establishment of combustion gas system and 
assumptions . 

2. In the temperature range from 200 to 1,500 K a modified 
form of the Bettie-Bridgman equation of state [3] is 
used to take into consideration the intermolecular 
forces at high pressures and lower temperatures. 

(The dew limits were not included in the expression.) 

3. In the temperature range from 1,500 to 6,000 K, the 
real gas properties were calculated from the chemical 
reaction equilibrium for individual thermal states 
(every 100 K in temperature, and at orders of 
magnitude of pressure from 0.001 to 1,000 kp/cm ) . 

Twenty different partial reaction products were considered 
in the combustion gas. 
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4. Using rational integral and transcendental polynomials, 
approximation equations were obtained for the enthalpy, 
entropy and the real factors, based on the derivations 
in (2) and the state points according to (3). 

5. i\ practical system of subprograms in FORTRAN is listed 
for application of the equations determined. With this 
system it is possible to compute changes in thermo- 
dynamic states immediately. 

Symbols Used 

Cp specific heat capacity at constant pressure 
F general thermod 3 mamic property 
H enthalpy (absolute) 

M molecular weight 
P absolute pressure 
R,Rq gas constant 
S entropy (absolute) 

T absolute temperature 
V volume 

Z real factor for the state equation Z = PV/RT 

n mole fraction of an individual gas in combustion mixture 

r air content in the combustion gas mixture 
i, j running indices 
y, empirical functions 
; air ratio, DIN 1940 

1. General Assumptions and Prerequisites 

The necessary number of independent parameters sufficient 
for adequate description of a system of combustion gases can 
be determined from Gibbs' phase rule: 

f = b - p + 2 (1) 

b = number of independent components in the system 
p = number of phases appearing 

f = number of degrees of freedom or the number of 
parameters required to describe the system. 
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In all combustion engines, the combustion gas appears 
only in the gas phase; that is, p = 1. The combustion gases 
under consideration here arise from oxidation of hydrocarbons 
in air» The air has a fixed composition except for humidity. 

Air humidity will not be considered here, though. Hydro- 
carbons of different compositions are used in combustion 
engines. In order to limit the range of computed points, it 
is necessary to select a fixed carbon/hydrogen ratio for the 
fuel. This ratio is taken here as C/H = 85.63/14,37 . This 
is exactly true for mono-olefins and naphthenes. From a 
comparison of various fuels for combustion engines, this 
average value seems suitable for a general representation. 

Such an assumed fuel, along with the air of fixed 
composition, makes up the two components of the total 
system; that is, b = 2. Then for the number of independent 
parameters, one finds that f = 3; that is, three independent 
parameters are required. 

For combustion gases, the following parameters are 
generally required as characteristic quantities and properties : 

a) Characterization of the type of gas 

Statement of the air /fuel mixture ratio; 
e. g. : 

air ratio . . , air volume before combustion 

stoichiometric air volume 

or 

F . . . Air /Fuel mass ratio 

volume of air in the combustion 

r . . . air content . . . gas after combustion 

total combustion gas volume 

b) Thermal parameters: 

Absolute temperature T; absolute pressure P and volume V. 
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c) Caloric parameters: 

Enthalpy H; entropy S; and specific heat capacity 
at constant pressure, Cp. 

Among the parameters in a), the air content, r, 
is especially suitable for an analytical representation, 
as Schiile suggested in the Z. VDI 60, p. 630, (1916) for 
fire gases. This is because the thermodynamic properties 
of combustion gases in the ideal gas state are linearly 
dependent on this parameter, due to the mixing rule. 

Aside from the generally common air content, r, the 
temperature and pressure must also be selected as independent 
variables. Then every combustion state with the assumed 
fuel is unambiguously characterized by the following 
parameters ; 

T absolute temperature 

P absolute pressure 

r air content 

Accordingly, the other parameters listed in a), b) and c) 
are dependent quantities. 

The relations between the dependent and the independent 
parameters are further represented by the thermal equation 
of state 




and the caloric equation of state 

H V-(T,P,r) 


( 2 ) 

(3) 


The thermal equation of sta.te represents the extra-molecular 
processes, and the caloric equation the intra-mo lecular 
processes. Both relations are empirical for real gases. 
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Because of the wide range of tne temperature and pressure 
for this representation, the chemical reactions at high 
temperature (dissociation) and the pres sure -dependence of 
the caloric properties must be taken into consideration by 
the thermal euqation of state at low temperatures and high 
pressures . 

If the empirical relations are known according to 
Equations (2) and (3), then the other thermod 3 niamic 
properties can be derived mathematically; for instance; 

S'- JdH.T or Cp- (dH/dT)p 


2, The Thermal Equation of State for Combustion Gases at 
High Pressures and Low Temperatures 

In order to consider the effects of real gases, 
euqations of state with several virial coefficients were used. 
The best- known generally applicable real equations of state 
are those of Beattie-Bridgman [3] and Benedict -Webb-Rub in. 

The Beattie-Bridgman equation of state is sufficient for 
the simpler technological gases. In its virial form it 
is as follows : 
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B(T), C(T) and D(T) are the second, third, and fourth 
virial coefficients. Z indicates the real factor, which 
has a value of Z = 1 for an ideal gas . Keyes , Smith and 
Gerry [4] changed the second virial coefficients of the 
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Beattie-Bridgman state equation in order to consider eqiaally 
the behavior of polar and nonpolar gases. The original 
Beattie-Bridgman equation is not suitable for representation 
of polar gases. The altered second virial coefficient of 
Keyes, Smith and Gerry is: 


B(T) Bo- 


AogDTJ 

Rf~ 


( 5 ) 


Stockmayer [ 5 j tested this form of the virial coefficient 
on binary gas mixtures: H2O-CO2, N2-NH2, and N2-H2O. The 

gases CO2 and N2 are not polar. H2O and NH^ are polar. 
Because of the polar water vapor content in the combustion 
gas, the altered second virial coefficient. Equation ( 5 ), 
should be used in the state equation. Equation ( 4 ). The 
constants for the pure gas components appearing here are 
known from measurements . 

The constants for the mixture are determined here 
according to the following equation 


in which the summations are to be carried out over all 
components of the mixture. 

n. , n. Volume proportions of pure gas components 

J- J 




Constants for the pure gas components 

Mixture constants due to interactions between 
two different types of gas molecules [ 6 ]. 


As yet there is no valid theory from which the 
constants K. . for i ^ i can be calculated from the 
constants for the pure components. But from various 
measurements and computations for binary mixtures [ 5 , 7 ] 
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bhe following relations can be used as approximations 
for the constants in the second virial coefficients of 
equation (5); 


A„ |A..-A„)'5 
B.. {J 8) (B. ” • B,/ 
{D. • D,,)’ •’ 


(7) 


Aitor application of these relations to combustion gas 
components in the range from pure air up to stoichiometric 
combustion gases, it could be established that the 
constants A, B and D in Equation (5) show an approximately 
linear dependence on the air content, r, in the form: 

K,„ = Kg + Kj^r. 

Because of insufficient reliability in the formation 
of the interaction constants K. . v^7ith 1 4 •] in the 
third and fourth virial coefficients, and because of the 
relatively minor contributions of the latter virial coefficients 
in the technological computations, only the first and second 
virial coefficients are considered for the combustion gas 
mixture. With the further approximation that W may be 
expressed as P/RT, one obtains the following initial 
solution for the thermal equation of state: 

V • Bir)- exp [(Dq • Du) T*) v 


Here the volume, V, is shown to be explicitly dependent 
on the independent parameters T, P, and r. 

With this equation of state, then it is also possible 
to obtain the pressure-dependence of the enthalpy, H, and 
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Uho cntj.*opy, S, by means of the following differential 
equations (b): 



Figure 1 » 


Specific l\eat capacity of the real stoichiometric 
combustion gas at pressures of 0,001 to 1000 kp/cm*^. 


I 

I 

KJ j 

piwi 






I--' 

Jj 

c 

w 


vW 




WX’' 


i,xm ‘ 

/ 

/ \ 


m 


ir«v' 


/ 

/ • ft ^ t 

i' ii 

i >J *® «V| 


t! 


a 


y f J 

/ - ! ,/■ 

^ ^ to 

■' , ; h» 

.'T iJ. ■ ' 


.\\V * W ' ^i\V * ta\V K 


Temperature T 


Figure 2. Enthalpy, H, of the real, stoichiometric 
combustion,. gas at pressures from 0.001 to 
1000 kp/cm-^. 


9 


ORIGINAL RAGE || 
OF POOR QUALITY 

3. Chemical Equilibrium Reactions in the Combustlo T ^ Gas 

Chemical reactions in the combustion gas must be taken 
into consideration in the temperature range above 1500 K. 

As the effect of the intermolecular forces decreases with 
rising temperature in this temperature range, and the real 
equation of state approaches the ideal one, the ideal 
thermal state equation can be used as the basis in this 
region. It is also apparent that the modified Beattie- 
Bridgman state equation selected also transforms into 
the ideal state equation at higher temperatures, even 
at relatively high pressures. At high temperatures the 
thermodynamic properties of a gas mixture are expressed by 
intra-mo lecular processes. The gas stores energy in the 
form of various states of motion, which lead to dissociation 
of the molecule into atoms or groups of atoms (radicals). 

An apparent equilibrium composition, dependent on T, P 
and r, is established between the reciprocal dissociating 
and bonding reactions. 

The calculation of the gas composition and, from that, 
the thermodynamic properties of the dissociating combustion 
gas are presented in [8] and [9] with respect to application 
of electronic data processing. The calculations performed 
here are closely based on [9] and are extensively described 
in [10]. The equilibrium composition must be determined for 
every thermal state point by means of a long iterative 
computation from the mass action law, the atomic balances 
(unchangeable number of the various types of atoms in 
different combintations) and Dalton's Law (sum of the 
partial pressures = 1) . 

In order to get a sufficient number of supporting 
values for the following approximations, the gas composition 
and all the thermodynamic properties of interest were com- 
puted every 100 K from 1500 to 6000 K at all orders of 
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magnitude of pressure between 0.001 and 1,000 kp/cra , and 

for the air contents (air ratios) r {;.)o(i); b,2(i,267);0,4 (i,7i3)! o,6 (2,603)j 
0.8 (5,278); 1,0 (x) 

This set of steps gave computations for about 2,000 state 
points . 

In the dissociating combustion gas, 20 different 
reaction products were considered as potential gas componcr.ts ; 

H 2 O, CO 2 , HCN, N 2 O, NO 2 , H 2 , C 2 , O 2 , N 2 , NO, CO, CN, 

OH, CH, NH, H, C, 0, N, and the noble gases of the air. 

The required absolute thermodynamic functions Cp/R, 

(H - Hq)/RT and S/R for the 20 individual gases in the 
ideal gas state, and the reaction enthalpies, were taken 
from [8] and [9], because they presented a critical selection 
of the latest data up to 1959. The data for the triatomic 
gases were from McBride and Gordon [11]. All the data used 
are tabulated in [ 10] . 

Of the results of computations from the reaction 
equilibria, only the heat capacity, Cp in kcal/kg K 
and the enthalpy, H, in kcal/kg, are shown here in 
Figure 1 and 2 for stoichiometric combustion gas (r = 0 
or A = 1 ). The absolute temperature in K is the abscissa 
for each figure. The seven curves presented in each figure 
correspond to the absolute pressures 0.001; 0.01; 0.1; 1; 

100; and 1000 kp/cm^. 

In these diagrams, two distinct reaction steps are 
recognizable from the wave-shaped course with rising temper- 
atures. These are due to the prevailing decompositions 
appearing in the various temperature ranges. The first 
step is for decomposition of the triatomic molecules into 
diatomic ones. The second step is for decomposition of 
diatomic molecules into monatoniic ones . 


4. Establishment of Approximation Equations 


okioinal 




In order to be able to work with the thermodynamic 
gas properties on programmed computer systems, these 
properties must be presented as mathematical equations. 

The following three functions are approximated: 

1. Real factor Z ^ >■) 

2. Enthalpy h-H(T, P, r) 

3. Entropy s s(t, p, r) 

The entropy is to be provided for the minimal system 
from the thermal state equation (2) and the caloric state 
equation (3) in order to avoid derivative errors in 
computing the entropy from differential calculus. 

Corresponding to the various effects of molecular 
processes on the thermod 3 mamic properties, one can 
decompose the three functions shown above into three 
parts each: ideal behavior, correction according to the 

real state equation (second virial coefficient) and 
dissociation correction. Accordingly, the three functions 
shown can be represented in the general form 


F •" Pldopl 'f* •! pRcalflos "P Foissojlolio 


( 10 ) 


Here the ideal gas properties, contain only 

the temperature dependence of the specific heat capacity in 
the undissociated ideal gas mixture. They are approximated 
by rational integral polynomials. The Tschebyscheff 
approximation procedure was used as the most practic£il 
method of approximation. 
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The correction term due to the selected second virial 
coefficients in the state equation, gas^ 

known analytically directly from the previous section, as 
are effects on the enthalpy and entropy from the differential 
equations (9). 

The effect of dissociation, AFj,. must be 

’ dissociation’ 

considered in the high- temperature range. The wave-form 
appearance of the dissociation contradicts an approximation 
by simple rational integral polynomials. For this reason 
it is advantageous to use transcendental types of functions 
for the approximation. For this case, a mathematical model 
was established by means of which the constants can be 
determined from a formula made up arbitrarily of trans- 
cendental and rational functions. On applying the principal 
of minimizing the squared errors, one gets with mixed 
transcendental-rational formulas non-linear determination 
equations for the constants. Tiiese can be solved by a 
Newtonian iteration. For this purpose, coarsely estimated 
numerical values are input into the computation for the 
constants which are initially unknown [10]. This approx- 
imation procedure makes it possible to match the most 
practical mathematical formulas to the curve form to be 
represented with enough free constants. 

After detailed studies and back-calculations with 
various types of functions for the curves of the real factors, 
enthalpy and entropy, the correction term ^F^igsociation 
was represented by the following formula; 

^ Foissoziaiion ai/(l+e'’' ) + a2/(l+e^2j (ll) 

In Equation (ll) the exponential functions y^^ and y 2 
and the amplitude functions a.-^ and a 2 can be solved with 
the aid of the curve to be approximated, specified point-wise; 
and they can further be approximated by the best-matching 
function types according to the method cited. 
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While and 72 are usable equally for the state equation, 
enthalpy and entropy, and a 2 must be matched separately 
to the particular thermodynamic properties. 




Applying the summation of Equation (10) to each of 
the thermodynamic functions to be represented here, one 
finally obtains as a result the approximation equations 
according to Table 1, which always contain the sequence 

^ “ '^^real gas ^ ^ideal ^ ^^dissociation’ real 

effect is superimposed in the ideal functions from the 

low- temperature side, and the dissociation effects from 

the high- temperature side. The approximation functions of 

Table 1 are shown in Figures 3, 4, and 5 for stoichiometric 

combustion gases. 


The three thermodyn-iamic functions are referred to 
the kmol of undissociated combustion gases, and made 
dimensionless with the ideal gas constant and the 

absolute temperature. By multiplication with the corresponding 
ideal gas constants, the thermodynamic data may be obtained 
immediately in the dimensions desired. The ideal molecular 
weight, Mq, of Table 1 must be used in relation to the 
weight (kg). 



Figure 3. Real factor for stoichiometric combustion gas 
from Equation (1). 
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Table 1: Approximation equations for the thermodynamic 

properties of combustioh gases. 


Table 1 


1 fieoKokior 
2(0in:rJ • • 

n 

^ £gJh.alpie : 


el ® B 




t HaSiiHHUL 
I « Ig'i 


, aociBlnn 




f/J 



a, -cviDscjt’t 

e, fO,m43» 

l•o,^7Jlt^ to,i$siU‘ r 

c, l•C,7l7i^l•lb‘OJHO^I•l'J■T 
ij ( • o^it3mi-o^oii)7')V' r 


Entro pig 


. JLf_d 


eL e 


^)\,t3,33ii-0,0fnrnnioooet Z<’,0' - Inn ♦ 

/.ff ' kto" 


-J,mr)n' 


i-" 


(UI) 


4 ml 


b„> >3,373 ns 

• 0,55J/<0 « f 

e/ • 0,BOS2I4 

y • V,5t)t37*lV 

* PAOSOOh r 

Cj» • 0,400 tOf 

V • 0,071 i49 

• o,s>‘fuj-;o' r 

eg* * 0,576 969 

b,> •0,4St50S 

- 0,3077,'0'll^ r 

Cj* • 0,060056 

0,109 164 

• o,i33ei.<id’- r 


^0,H9 75i>l0 



*0,600:J4^IO 

• 0,333073 W‘ r 





• 0,3060 

* to, 573 • 0/303 In 


so?«-^j . n "(I' r/(j;;w}o,Mss'if‘’^' * 


m 10,0371 - J,0 e ii'*'’" 

A !T 3,77305'ib*~ o,90O7it>w‘ r , 0,. ff,s««7no' r ; 0 * OjOOOOoy- qaohii' r 

5 Temperolur !0 A DfUCk 

g ..T n , .TF.lRteiy) • n . P .. . P olm 

1000 K lOOO'(yg) «R ' 1 ato 0,9B066B bcr 0,967641 atm 

7 Luftoeholt r: ■ it Q 8 G°s*<°nst°n<c: R „ 

' \ — r ; Oaskonstqnle des Idaolen Gases X 1 

1,069B.X-1 VGu„dj„„i,.ri 2^ Q 

X 2 Molmossed. ideolen VG =28.8975S*0, 06021 r[kg/kmol), 


1 Real factor 

2 Enthalpy 

3 Entropy 

4 with 

5 Temperature 

6 Pressure 

7 Air content 

8 Gas constant 

9 Air 

10 Combustion gas^„j. 

11 Gas constant for the ideal gas 

12 Molecular weight of the ideal combustion gas: 


15 


Figure 4 



Reduced enthalpy for stoichiometric combustion 
gas according to Equation (2). 



Reduced entropy for stoichiometric combustion 
eas according to Equation (3). 


Figure 5. 



Figure 6 
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Error curves for enthalpy. 



Figure 7. 


Error curves for entropy and realf actor , 
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I 
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common ABi AC iGiJiVliLPiPIlfCfM 
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COMMON VIOi V 20 i V 1 i IS 

NA « 0 

CAU ANF 

UTA • T • ,001 

CAU. PPI 

IF INAI litOil 

10 CALL VA 

If INAI Ifllil 
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CNO 


SUBMQUTINC HATP 
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CALL PPI 
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END 


SUBROUTINE SATP 
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COMMON Ti P| Vi Hi Si Mi NA 
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COMMON BZIM«e3n)|Cl8liK3IA)i03|a| 
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NA » 0 
CALL ANF 
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IF INAI liUil 
12 CALL CMISCM 
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END 
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EXTERNAL VA 
NA r 0 
G • V 

CALL T/UMIVA) 

return 

ENC 


SUBROUMNF TAPH 
REAL J« Kl, K?i K3i LPi MS 
COMMON T, P| Vi Hi Si Ri NA 
common ABiACiG.JiVltLPiPItTriA 

CCKMON B2l7he3miCia)iK3IAI|03(fl) 
COMMON VlOi Y20f Y3| YA 

external ha 

NA • 0 

G • H 

call T/UHIHAI 

RETURN 

END 


SUBROUTINE TAPS 

RFAL Ji Mil K 2 » K 3 i LP» M$ 

COMMON Ti P| V. Hi Si R« NA 
COMMON AEi ACiGi JiVliLPiPli TEIA 
COMMON fl2(7lie3(7liC| eiiK3(AI|Q3(fil 
COMMON YlOi Y2Q, Ti. VA 
EXTERNAL SA 
NA •- 0 
G - S 

CALL T2UMISAI 

RETURN 

ENC 


SUBROUMNC PATV 

REAL Ji Kl. K2, K2, LP| MS 

COMMON Ti Pi V. H. Si H« NA 

COMMON AB.ACiCt J.VliLPiPli TETA 

COMMON B2l7lini7liCl BliKH A) ,03161 

COMMON YIC, Y2Q, VI. YA 

external VA 


subprogram system. 


SUBROUTINE PATS 
REAL Jf Kli R2f K3i lPi MS 
COMMON ft Pi V. Hi Si Ri NA 
CORPON A«|ACiGiJfVliLPiriiTrTA 
COMMON B2ITI|I3I?I«C(B||K3|AI,03IBI 
common YlQi Y2Q| V)i YA 
EXTERNAL SA 
NA • 0 

IF ISI li It 3 
1 G > ALOCISI 
CAU PIUMISAI 
IF INAI 1|2|2 

1 NA « I 

2 RETURN ^ ) 

INC 


SUBROUTINE fPAVH 
REAL Ji Kit K2i M3i LP| MS 
COMMON Ti P| V| Mi S. R. NA 
COMMON ABiAOiGiJiVTiLPiPI iTF TA 
COMMON B2l7li B1l7}tC(BliM3|A)rQ3TBt 
COMMON VtOi V?Oi V3» YA 
external ha 
NA - 0 
J ■ H 

CALL TPZUMIHAI 

RfTURN 

CNO 


subroutine tpavs 
real Ji Mil M2i K3« LP« MS 
COMMON Ti P| V| M. $. R| NA 
COMMON ABiADiCt Jf VtiLPfPt iTFTA 
COMMON B2<7lf B3l7l|C(ai|K3IATi03(6l 
COMMON YlOi Y20| Yli VA 
EXTERNAL SA 
NA • 0 
J * S 

CALL TPZUHISAl 

RETURN 

ENC 


FUNCTION UFTCCHICI 
COMMON T| Pi Vi Mi S| Ri NA 
NA • .0 

IF <E~UI Ii2i2 

2 IF EE - lOOOU 3iAiA 

3 R • lE-t.l/l .C69B*E) 
tFTCEM ■ R 

CO 70 A 
A R*U 
IFTGEH-R 
(0 10 5 
1 NA • I 
5 RETURN 
CNO 


SUBROUT INE VA 

REAL Ji Kli K2i K7i LP| MS 
COMMON Ti Pi Vi Hi Si Ri NA 
COMMON A6, ACfOi J|VT|LP(Pli TETA 
CCHHCN R2l7),eil7)iC(ei.K1EAIrOT(ai 
COMMON YlOi Y20, Y3i VA 
RO * B3IA.U 
CALL CR 
IF (NA| t.lO.l 

to A » OM I !•?!*•. OOTi'TlCUOlE?) *0)1)1 
IHPI/VaCtU ♦ IK1I 2)-tAJ»PI /TCIA 
Xn - A*R0»T/EP*K3IAT| 

1 RETURN 
(NC 


SUGROUTINE HA 

REAL Ji KU M2| K3i LPi MS 

COMMON Tt Pi Vi Mi Si Ri NA 

COMMON AeiAO|G|J|Vl|LP«Pt|TnA 

COMMON R2(7)ie3l7),C|6>|K3{A),Q3(B) 

COMMON YlOi Y2Qi Y3i YA 

RO « BIlA.ie 

CALL GP 

IF INAI lilCfl 

10 f • HRNIfl2.rETA,7J*l03(AT*PI**,0n5/ 
IYIO*IOHSI-3.04A*LPI^Y20*|K?U»-?, 
2«VA*| l.«Y3l MPI UTCTA 
AO - RO»T*«/K)(A» 

I RETURN 
CNC 


SUBROUTINE $A 
REAL Ji KIi K?i K3| LPi MS 
COMMON Tf Pi Vi Mi Si R. NA 
COMMON A8« ACiCi JiVliLPiPI iTETA 
CCKMON B2( 71i B3m,CI BiiMI AI.O3I0I 
COMMON Y!0| T?0, Y3, YA 
PO » B3IA.18 
CALL CR 
IF (NA) liTit 
3 TF IT) 2i 2, A 
A IF (PII2, 2. ‘ 

S AC » 101l6>*ALnC|TI*HRNlB3.Tb*TA«7)* 
UC1l7)»Pl**l-.01A52l)/YlQ-LP*l03lfl) 
2nA.5l?*PI**l-.OSfcAB) I/Y20-IYAAPI* 
1U.*Y3I )/TUA)*R0/K3IA» 

GO TO I 
? NA » I 
I RETURN 

END • - ,, J 
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SUBROUlINt ANF 

REAL Jf Kii K2> Kli LPi MS 

COMMON 1 1 Pi Vi Hi Si Ri NA 

CCMPQN ABiAO|CiJiVl|LP|PI|TE?A 

COMMON R 2 miC 1 (T)iC(BtiR 1 IAIiQ)fBf 

COMMON VIQi V 7 Qf Vli VA 

ClMfNStON KlIAIi 71A), 01<BI| 02161 

Rim ■ 2.77icsr-» 

KH2I - 6.A22171 
RU3) « 0.4B6IE- 
KUA) • 2e«a979B 
K2IU • -iSCOTlU-A 
R2I2I • -.983ATC-A 
K2I3I • -.Al3lf-2 
K2IA) • 0.06021 
CALL KOMBlAi Kli K2) 

CO 10 I - li A 


10 

F)l II 


cm 



eti II 


0.A2C 



(1121 


0.77A 



CIOI 


-0.012a 



CM«I 


25.9 



CUSI 


66. S 



CUfcl 


3.SS66 



cim 


8.179 

'r 


cum 


0. 

V 


(2111 


-0.193 

1 


(2121 


-0.U9 



C2I]I 


-o.oos 



(2141 


-ll.O 



(2151 


-IA.35 



caut 


-0.0659 



C2ITI 

« 

-3.72E 



(2181 

* 

-2iA3A 



c«u 

KCH8U, «l> (21 



cn M 

1 « ii a 


11 

cmi 

• 

cm 



AB • -O.9C06 • I-.2T72 • Vl*«A * VI) 

RETURN 

(NC 


SUBROUTINE KCMBlNi A, 0) 

REAL Ji Kli K2| K3i IP . MS 
COMMON Tf Pi Vi Hi S( i. NA 
COMMON ABiACiCiJiVliLS^iPliTETA 
common e2n)tB3mfCia)iK3(A)|Q3|6) 
common YlOf Y2Q, V3i VA 
CtPFNSION Aieii BUI 
CO 10 I • li N 
10 cm • Ad) » BUMR 
RfTURN 
ENC 


subroutine PPI 

REAL Ji Kit K2| K3i LP» MS 
COMMON Ti Pi V. H. S. Ri NA 
common AeiAOiCiJiVliLPiPliTETA 
COMMON R2mi e3l7)|CiaiiK3IATi01(ai 
CCMMOM YIO* Y2Q, Y3, YA 
IF (PI 111. 3 
3 PI - PM.C1S72E-5 
LP « ALOGIPl ) 

CC TO 2 

1 NA » 1 

2 RETURN 
f»P 


SUBROUTINE CR 
REAL Ji KU K2. K). IP. MS 
COMMON Ti P. V. H. S| R. NA 
COMMON ABi APiGiJiVliLP.PI iTETA 
COMMON G2(7Tie3(7l*Ciei,K3IA),Q3(81 
COMMON YlCi Y20i Y). VA 
IF (TETA) 1, I, IC 
10 f * ALOCI 1ETA/U6ST 

A ■ -.5l*PI*d-.lC3)-.lA ♦ 0«(.l2- 
U2S«PI«*|-.127)) 

A « l5.«tXP(A«B)i(.573i.oe3HP)» 
Uxp|-OC6«TETA)-6,907e 
Y| » AB«Pl»«(*.02l?)*A 
Y2 • Ift. 0972-2. A3*IETA*PI**I-.0AS) 
YIC • l.iEXPOn/TETA 
Y2C » l.iEXPn?I/TCTA 
Y1 » K3m/ITETA*!ETA) 

YA ■ EXPIY3JAK3I ll/TEIA 
CC TO 2 

1 NA . 1 

2 RETURN 
ENC 


SUeRCUTINF CMlSCH 
PEAL Ji Kit K2« K). LPi MS 
COMMON 1 I Pi V« Hi S. Ri NA 
COMMON AeiADiGiJ.VlfLPiPIfTriA 
COMMON R2l7)ie3<7liCIB)iK3l A)|Q3(B) 

COMMON YlOi Y2Q, Y3i YA 

CIMfNSinN Midi, H2I7I. SU7)i S2(7I| MSIA) 

Him - 3.S1A9S6 

H)(2) « 0.131A22 

Kl<3) » O.A7m'2 

HUAI - -.28736T 

HUS) • 0.7A256 1E-T i 

HU6I • -.9U3AAE-2 
HU7I - 0.A39e9Ef-‘1 
Ham • -.Q0SC2E 

K2I2I ■ -.38350A » 
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1121 11 

» c.ieiiiA 


H2IAI 

» -.t9A8e.'f-i 


H2I5I 

> O.UACAIt’l 


H2I6I 

• *.2CA53IE-2 


H2I7I 

• O.ICIMCE-3 


CiiU 

KCMBI7, HI, H21 


CC IC 

1 * 1, 7 

10 

82III 

• Cl 1 1 


smi 

• 2.912979 


StUi 

• 0.583BJ7E-1 


SU3I 

• 0.87I3AS 


sue 

■ -.A5I556 


SM5I 

• 0.1C9UA 


SH6I 

. -.12975AE-1 


S1I7I 

• 0.6C823AE-] 


S2I 11 

• 0.5S3IAC 


S2I2I 

- -.AB5CCI 


52131 

• 0.698ielE-l 


S2IAI 

. -.287720E-2 


52151 



S2I6I 

« C.9202’ir-A 


52171 

• -.299£79E*5 


CALL 

KCM8I7. SI, S2) 


CO 11 

1 > 1, 7 

II 

£311) 

• cm 


MSI II 



MSI?) 

• -.ACC98I 


M5I3I 

• 0.576989 


MS(A| 

. C.C60056 


8311) 

« a?llllMMN(M5, 


TURN 

ENC 


FUNCTION HPNIBI, X, 



Jf Klf Kjf 


COMNDN r, Pi V, Ml Si Pi NA 
COMMON PBiACiCiJiVliLPiPliItIA 

common e2i7i,e3miCie)iAAiAiiCiiai 
COMMON yiOi yjo, vji 
CIMENSION 111171 
HPN • Oi 
I • N»l 
I ! » 1-1 

HPN r HPN • A • BUM 
If 11-11 Jl ?i 1 
? PcroPN 
ENC 


SIBPCUTINE I20MI1KI 
PEAt Ji Kli A7i M’l LPi MS 
COMMON I, Pi V| Hi Si Pi NA 
COMMON ABiACiGiJiVliLPiPliIEIA 
COMMON e7l71iE?(71iCieiiM3IAIi03l81 
COMMON »lOi y?Qi Y3f YA 
EKtEPNat IP 
CALL CMISCH 
CALL ANF 
CALL PPI 
IF INAI Ii17iI 
12 T « OLUiaAPiiAOOO.ilPiIM 
NA « ABSINAI 
1 PE.IUPN 
ENO 


StDPDUIlNF PEUMIPKI 

REAL Ji Ml, K?| K2i LPi MS 

CCHPCN Ii P| Vi Hi Si Pi NA 

COMMON AUiAOiCiJ,VliLPiPI|IFIA 

COMMON D2i7iie3i niCieiiK3iAiiC3iei 

COMMON YlCi Y70i Y3i YA 

EYICPNAL CP 

PC • asiA.iB 

CALL CMISCH 

CALL ANF 

t£IA . I • .001 

IF (I-AOC.I *iAi3 

A YA-E>P|A3m/IUIAMEIAI|.K1l ll/IFIA 
P • RO*1/IIV-IK3l2l-YAl»aA. 78308/ 
IK3IAII»(3IAII 
IF IPI SfSi6 
6 PI ■ PP1.01S77C-5 
LP « ALOCIPII 
CO 10 I 
5 NA * 1 
CO 10 I 

3 LP - ULLin.ilAiHiDPiPM-lO. 

IF INAI li7il 

2 PI • EXPILPI I 

P » P|*.Ve066SE5 
1 return 
ENO 


SLePOUTINE IP2UNIW1 
REAL J, Kli K2i K’l LPi MS 
COMMON T, Pi V, Hi Si R| NA 
COMMON AeiAOiGiJiVliLPiPliTEIA 
COMMON B2l7lie3l7l,CieliK3IA|.03l8l 
COMMON YlOi Y20i Y’l YA 
EAIEPNAL FIF 
IF IVI 1. Ii 3 
3 C- • ALOCIVI 
CALL ANF 
CALL CMISCH 

T ■ ULLII2AS.i7CC..FIFiM 
IF INAI Ai2i2 
A I « ULLU 7CC.i20CC.iFIF.lil 

IF INAI ;.2i2 I 

5 T ■ ULLII 2COC.i6CCO.iFIFiVI 1 


1 

IF INAI ti2.; 

1 NA • I 

2 PtlLPN 
INC 


ILNCItON CPIAt El 
PEAL A. All M?I LPi MS 
COMMON r. P. V. Ml $. Mr NA 
COMMON Ae,A0iCi3iVi.LPiPlilE|A 

COMMON BaniiEsi n.ciB).M 3 iAi.e)iei 

COMMON YIB. Y?0| Yli YA 
LP » «-IO. 

PI • EXPILPI 
F • PI • isecetSES 
CALL E 

IF INAI 2.3.2 
3 IF lACI I. 1. A 
A CP « alOGIACI - C 
CO 10 2 

1 NA ■ 1 

2 RtIURN 
INC 


FUNCTION TPIX. El 

REAL Ji Kl. K2. KJ, LP. MS 

COMMON t. P, V. H, S. R. NA 

COMMON AP|AC|C.J,Vl,LP|PI|irtA 

COMMON 82171, B3I7I,C1PI,K3IAI,C1IBI 

COMMON YlCi Y70, >3. YA 

I • X 

UIA • I • .001 
CALL F 

IF INAI 1,2,1 
2 IP . AO - 0 
1 RFIURN 
INC 


FUNCIION FIFIX.EI 
REAL J, Kl, K 2 , K 3 , LP, MS 
common t, P, V, II, S, P, NA 
COMMON Ae,AC,G,J,Vl,LPiPIilEIA 
COMMON E 2 ( 7 l,E 5 l 7 l,Ciei,K 3 IAI, 03 iei 
COMMON YIO, Y 25 , Y 3 , YA 
EXIEPNAL OP, VA 
PC • B 31 A. 3 E 
I » X 

lETA • I • .OCl 
IF ll-ACC.I SiS,l 

5 YA«CXPIK 3 I 3 I/I 7 EIA»IEIAI |*K 3 I il/IEIA 
P » PO*l/HV-lK 3 l 2 l-YA)*eA. 78308 / 

IK 3 IAIIPK 3 IAII 
If IPI 3 ,!it 

6 PI . P*l.CIG 72 E -5 
LP « ALOCIPII 

CO TO 7 

1 ARC « RC»T/IVA.SBOttSE 5 *K 3 IAI I 
IF lARCI 3 , 3,2 

2 AA • A.SStA’.CCIAPGI 
EP. ■ AAI.EA 

LP » ULL 2 IAA, 8 B,DP,YAI- 10 . 

If INAI A,E,A 

8 PI - EXPILPI 

P • Pl». 980 etif 5 

7 CALL E 

IF INAI A, 9 , A 

9 FIF - AO-J 
CO TO A 

3 NA ■ I 
A PFIUPN 

INC 


FUNCIION ULLIlAi 8, F, El 

BEAL J. Kl, K2, K3, LP, MS 

COMMON I, P, V, H, S, R, NA 

COMMON Ae,AC,C,J.v),LP,PI,1E1A 

COMMON f2l7lie3ni,CI8l,K3IAI,C3iel 

COMMON YK, Y20, Y3, YA 

NA • 0 

N2 • 200 

EPS • 0.CC5 

XON ■ A 

XAN ' 0. 

YON » FI XON, El 
IF INAI lAilSilA 


15 

X2N 

« B 



Y26 

• FU2N 

. El 


ir 

(NA) 

29 

ir 

(Y2NI 

27,25, 2e 

2b 

ir 

1 YON) 

12,28,30 

27 

IF 

(YONI 

30,28, 12 

12 

» 

N2 - 1 



IF 

INZ 1 1 

3,13,16 

16 

XIK 

• IXON 

• X2NIP.5 


tlN 

■ FIXIN 

, El 


IF 

IHAI !A,17iU 

17 

DN » 

Y2N - 

YIN 


IN • YIN-YON 
CN • ON-EN 
ON • ICN«ENI«.S 
ON ■ 2.»YIN*CN 
IF ICNI 18,1,18 
16 ARC • CN»ON - ON 
IF lARGI 13,19,19 

19 HN • SORTIARCI 
EN - IHN-DNI/CN 

IF lAeSIENI-l.n 2,2,20 

20 EN - -ICN«HNI/CN 

qc TO 2 


1 IN » -Y IN/ON 

2 XIN • ENP|X1H-I0N|*X1N 
YIN » FIX3N, El 

IF INAI IA. 2 I.IA 

21 IF lYINI 3 i 22 , 22 

22 IF IY 2 NI A, 4 , i 

3 IF IY 2 NI t, A. A 

S Y 2 N • YIN 

X 2 N • XIN 
EC fC 6 
A YON . yin 
XON . XIN 

8 IF IY 3 NI T, 23 , 23 

21 IF lYONI 8 , 9 , 9 

7 IFIYONI 9 , (. 8 

9 YON • Y 3 N 
XCN • X 3 N 
CO 10 10 

a Y 2 N • YIN 
X2N • X!N 

10 IF IABSIX 3 N-XANI-EP 5 I It. 2 A, 2 A 
2 A XAN • X 3 N 

CO 1 C 12 

11 ULII • X 3 N 
CO TO lA 

25 ULLl • X 2 N 
CO 10 lA 

26 ULLI • XON 
CO 10 lA 

30 NA • -I 
CO TO lA 
13 NA « I 
lA REIURN 
IND 


29 IF 

26 If 

27 IF 


FUNCTION ULL2IA, B, F. El 

REAL J, Kl, K2, K3, IP. MS 

COMMON T. P. V, H, 5, R, NA 

COMMON AEiACiC.J.YI.LPiPI iTETA 

common B2I7I, 83171, Ciai,K3IAI. 03181 

COMMON YIO, Y2Q, Y3, YA 

NA • 0 

hi « 200 

FPS » .005 

XON • A 

XAN • 0. 

YON ■ FIXON, El 
IF INAI 1A,1S,1A 

15 X2N ■ 8 

Y2N ■ FIX2N. El 
IF INAI tA.29.IA 
IY2NI 27.25.26 
lYONI 12,28,30 
lYONI 30,28,12 
12 N2 ■ Ni - 1 

IF IN2I 13,13.16 

16 XIN • (XON*X2N)»,5 
YIN » FIXIN, El 

IF INAI lA.n.lA 

17 CN • Y2N-Y1N 
EN » YIN-YON 
CN ■ CN-EN 

ON • IDNTENM.5 
CN • 2.*Y1N*CN 
IF ICNI 18, I, 18 

18 ARC • ONACN-GN 

IF lARCI 13, IS, 19 

19 FN - SORTIARCI 
EN « IMN-CNI/CN 

IF lABSIENI-l.ll 2. 2. 2C 

20 EN • -IDNAHNI/CN 
CO TO 2 

1 EN • -YIN/ON 

2 X3N ■ EN*IX1N-X0NI«XIN 
Y3N » FIX3N, El 

INAI IA,21,IA 


IF 

21 IF 

22 IF 


lYlNI 
IY2NI 
IF IY2NI 

5 Y2N . YIN 
X2N • XIN 
CO TO 6 

A YON - YIN 
XON • XIN 

6 IF IY3NI 
23 IF lYONI 

7 IF lYONI 

9 YON . Y3H 
XCN « X3N 
CO 10 10 

8 Y2N • Y’N 
X2N • X3N 


A, 


7 

8 

9, 8 


22 , 22 


23, 23 
9, 9 


IC 

IF 

lABSl x; 

29 

X9N 

' ■ 

XTN 


CO 

TO 

12 

11 

ULL2 • 

X3N 


CO 

10 

19 

25 

UU2 - 

X2N 


CC 

TO 

19 

2fl 

UU2 - 

XON 


CC 

TO 

19 

10 

NA 

■ - 

1 


CC 

TO 

19 

13 

NA 

- 1 


19 

RETURN 



II, 2A, 2A 
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Table 2, Organization of the THERMOBIBL subprogram system. 


Serial Subprogram Independent Dependent Approximation Equation 
No. name variable variable “llblf 1® 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 


VATP 

HATP 

SATP 

TAPV 

TAPH 

TAPS 

PATV 

PATS 

TPAVH 

TPAVS 


V 

H 

S 

T 

T 

T 

P 

P 

T and P 

T and P 


UFTGEH U) 


from: ; 
T,P,r 

T,P,r 

T,P,r 

P,V,T 

P,H,r 

P.S,r 

T,V,r 

T,S,r 


V.H.r 

V',S,r 


air X ratio 


indirectly from I 
indirectly from II 
indirectly from III 
iteratively from I 
iteratively from II 
iteratively from III 
iteratively from I 
iteratively from III 
iteratively from I and 
iteratively from I and 


II 

III 


By comparing the values from the approximation equations 
with the base data previously obtained pointwise, one 
gets the errors of the approximation. These are shown 
in Figures 6 and 7 for the real factor, the enthalpy, and 
the entropy. The errors are greater than in the usual 
approximations of simpler functions, corresponding to the 
complexity of the functions to be represented and the 
difficulty of the three-fold dependence of T, P and r. 

In selecting types of approximation functions and the number 
of free constants in them, one must make a compromise 
between the accuracy of the approximation and the computing 
cost and storage capacity of the computer. From the 
diffuse course of the error curves, it is also apparent 
that the equations established represent a certain optimum 
between smoothness and the exact reproduction of the already 
scattered basis data determined point-wise by iteration. 
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Comparison of the thermodynamic data computed from 
the equations of Table 1 with those from the Mollier I,S 
diagrams of Pflaum [2] shows only errors of less than 2% 
in the entire range presented by [2], The diagrams [2] 
themselves have been compared with other modern thermal 
diagrams in [12] with very good results. 

5. Application of the Approximation Equations in Programmed 
Computer Systems 

By means of the approximation equations, Table 1, 
a system of subprograms was established (THERMOBIBL) . It 
provides for computing the explicit variables V, H, and S 
from T, P and r, as well as computation of the inverse 
functions; e. g., T from P, S and r. Solution of the 
inverse functions is possible only iteratively. In these 
subprograms, one of the variables T, P, V, H, or S can be 
calculated from any three of the others, as the scheme of 
Table 2 shows. The air content, r, is always used as an 
independent parameter, corresponding to the general usage 
of the I,S diagram. 

"Quadratic bracketing" according to Zurmuhl [13] 
was used as the iteration procedure to solve the inverse 
functions. The Newtonian iteration procedure failed because 
of the wavy form of the approximation functions, Figures 
3, 4, and 5. 

The subprogram system which provides for computing 
the dependences according to Table 2 with call of the 
particular subprogram desired is summarized by the name 
"THERMOBIBL" and was written in FORTRAN (DOS basic FORTRAN IV). 
It has been listed in Table 2 for immediate use. This program 
system is a further development of that published in [14]. 
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All the variables appearing here have units from the 
international MKSA system: T(K); P(N/m^)j V(m^/kg); H(J/kg); 

S(J/kg K)j r (dimensionless). The ideal gas constant used 
is : 

Rq = 8314.38/28.898 + 0.06 r) in J/kg K 

With respect to beliavlor of the equations at low 
temperatures and high pressures, the state range, which 
is represented by the iteratively obtained inverse functions 
(Table 2, No. 4 - 10) is limited to the following boundaries: 


T - 250 - 6,000 K ^ 

P = 0.001 - 300 kp/cm" 


98.0665 - 29.4-10^ N/m^ 


There is an additional pressure limit for TPAVH and TPAVS: 

2 

Maximum pressure 200 kp/cm in the range 250 - 400 K. 


The following details are mentioned for setting up a 

computer program using the THERMOBIBL subprogram system. 

1. The following COMMON region must be specified in the 
main program: 

REAL J, Kl, K2, K3, LP, MS 

COMMON T, P, V, H, S, R, NA 

COMMON AB, AD, G, J, VI, LP, PI, TETA 

COMMON B2(7), B3(7), C(8), K3(4), Q3(8) 

COMMON YIQ, Y2Q, Y3, Y4 

Here only the first seven quantities of the COMMON 

are outside the subprogram of interest. 

2. The variables established in the COMMON region have 
the dimensions specified above and are REAL quantities 
except for the INTEGER, quantity NA. 

3. The entire subprograni system, THERMOBIBL, should always 
be stored in the computer. 
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4. There should be no new definition of the following 
names in the main program, because they appear as 
program names in the subprogram system: aside from 

those listed in Table 2, the names VA, HA, SA, ANF, 
KOMB, PPI, GR, CMISCH, HRN, TZUM, PZUM, DP, TP, FIF, 
ULLl, ULL2. 


5. Before the call of one of the subprograms listed in 

Table 2 (e. g. , CALL TPAVH) the input parameters 

(independent parameters, e. g., V, H, r) must be 
placed in the COMMON region. After running the 
subprogram the corresponding dependent variable 

(e, g. , T and P) must be fetched from the COMMON 
region. 

6. After running one of the subprograms listed in Table 2 
it is possible to test for the correct course of the 
computation via the quantity NA: 

a) NA = 0 correct course of computation 

b) NA « 1 computation erroneous because the range 

of presentation has been exceeded 
By means of this control quantity it is possible to 
branch or terminate the main program corresponding to 
certain requirements without destroyiiig the computation. 

This subprogram system makes it possible to compute 
all the thermodynamic quantities with combustion gases rapidly 
and accurately on a computer. The simple and general appli- 
cability of the subprograms makes research and design work 
easier for combustion machines. 


ORIGINAL PAGE IS 
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